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ABSTRACT 


Recent S'piteer/InfraRed Array Camera (IRAC) photometric observations have re¬ 
vealed that rest-frame optical emission lines contribute signficantly to the broadband 
fluxes of high-redshift galaxies. Specifically, in the narrow redshift range z ^ 5.1 - 
5.4 the [3.6] — [4.5] color is expected to be very red, due to contamination of the 4.5 
fim band by the dominant Ha line, while the 3.6 /im filter is free of nebular emission 
lines. We take advantage of new reductions of deep Spitzer/IRAC imaging over the 
GOODS-North-bSouth fields (Labbe et al. 20151 to obtain a clean measurement of the 
mean Ha equivalent width from the [3.6] — [4.5] color in the redshift range z = 5.1 - 
5.4. The selected sources either have measured spectroscopic redshifts (13 sources) or 
lie very confidently in the redshift range z = 5.1-5.4 based on the photometric redshift 
likelihood intervals (11 sources). Our Zphot = 5.1-5.4 sample and Zspec = 5.10-5.40 
spectroscopic sample have a mean [3.6] — [4.5] color of 0.31 ± 0.05 mag and 0.35 ± 
0.07 mag, implying a rest-frame equivalent width EW(Ha-|-[NH]-b[SH]) of 665 ± 53 


A and 707 ± 74 A, respectively, for sources in these samples. These values are con¬ 
sistent albeit slightly higher than derived by Stark et al. 2013 at z ~ 4, suggesting 
an evolution to higher values of the Ha-|-[NH]-b[SH] EW at z > 2. Using the 3.6 pm 
band, which is free of emission line contamination, we perform robust SED fitting 
and find a median specific star formation rate of sSFR = I 7 I 5 Gyr“^, 7 ^ 2 ^ higher 
than at z ^ 2. We find no strong correlation (<2 (t) between the Ha-|-[NH]-|-[SH] EW 
and the stellar mass of sources. Before the advent of JWST, improvements in these 
results will come through an expansion of current spectroscopic samples and deeper 
Spitzer/IRAC measurements. 


Key words: galaxies: high-redshift - galaxies: formation - galaxies: evolution 


1 INTRODUCTION 

In recent years, large multi-wavelength photometric surveys 
conducted with the Hubble and Spitzer Space Telescopes 
have enabled us to study the properties of galaxies over cos¬ 
mic time. Synthetic stellar population modeling of broad¬ 
band spectral energy distributions (SEDs) has led to the de¬ 
termination of various physical properties (e.g., stellar mass, 
star formation rate, age, dust extinction) of these galaxies 


(Eyles et al.|2005 

Yan et al.|2006 

Stark et al.|2009 Gonza- 

lez et al.|2010[ IBouwens et al.|2009 20121 Finkelstein et al. 

2012; Oesch et al.|20131 de Barros 

et al.|2014| ISalmon et aL| 


2015). Of all derived quantities, stellar masses are particu- 


age of the stellar populations, metallicity or other parame¬ 
ters have no significant effect on the estimated masses (e.g., 
Finlator et al.|2007 Yabe et al.||2009 ). 


Considerable progress has been made in refining current 
characterization of galaxies from the observations. Even so, 
the measurements of the specific star formation rate (sSFR, 
the star formation rate divided by the stellar mass) have 
presented a puzzle to the theoretical understanding of the 
build-up of mass in galaxies, (e.g., Bouche et al.|20i0 Wein- 
mann et al. 2011). Several past studies had indicated that 


larly robust in stellar population fits. Small changes in the 


the sSFR of sources with fixed stellar mass shows no ev¬ 
idence for significant evolution between z ~ 2 and z ~ 7 
(Stark et al.|2009 Gonzalez et al.|2010). This result was in 
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apparent disagreement with semi-analytic models predicting 
a strong increase in the specific inflow rate (i.e. inflow rate 


divided by halo mass) of baryons with redshift (Neistein & 
Dekel|2008| ). 


Subsequent work has strongly suggested that the sSFR 
of galaxies is somewhat higher at 2 > 2 than it is at 2 ~ 2 
(e.g., Schaerer & de Barros 2010; Stark et al. 2013; Gon¬ 
zalez et al. 2014; Salmon et al. 2015). The change in the 
inferred sSFR evolution with cosmic time was the result of 
improved observational constraints and a more sophisticated 
treatment of those constraints. One example of this is in a 
consideration of dust extinction in computing the sSFRs at 
2 > 2. While it was not possible to properly account for 


the impact of dust extinction in initial work (e.g., Gonzalez 
et al.|[2010l due to large uncertainties on the UV colors of 


2 > 2 galaxies, later work (Bouwens et al. (20121 was able to 


make use of new measurements of the UV continuum slopes 
to account for the impact of dust extinction, finding ~ 2x 
larger sSFRs at 2 ^ 4. This suggested a 2x evolution rela¬ 
tive to the 2 ~ 2 value ( Reddy et al.|20l2 |, but still leaving 
the sSFR at 2 > 4 approximately constant. 

Even more important has been the increasing aware¬ 
ness of the impact of rest-frame optical nebular emission 
lines (e.g., Ha, [OIII], [Oil]) on the broadband fluxes (e.g.. 


Schaerer &: de Barros|[^09[ [2010 Schenker et al.|[2013 1. At 
high redshifts these emission lines are shifted into the in¬ 
frared, contaminating the IRAG measurements of the stel¬ 
lar continuum. Inferred stellar masses from fitting of stellar 
population models will then be overestimated, resulting in 
an underestimate of the sSFR. 

Since these strong rest-frame optical lines are inacces¬ 
sible to spectroscopy beyond 2 ~ 2-3, the strength of the 
nebular emission lines has been estimated from the contam¬ 
ination of the Spitzer/IRAC 3.6pm and 4.5pm bands for 


galaxies at 2 > 3. Shim et al. (20111 show that galaxies in 


the range 3.8 < 2 < 5.0 are considerably brighter at 3.6 pm 
than expected from the stellar continuum alone and argue 
that this excess is due to strong Ha line emission. [Stark et| 
al. (20131 derive an Ha equivalent width (EW) distribution 
by comparing the [3.6] — [4.5] color of spectroscopically con¬ 
firmed galaxies in the redshift range 3.8 < 2 < 5.0, where 
the Ha line lies in the 3.6 pm band, with an uncontaminated 
control sample at 3.1 < 2 < 3.6. The results indicate a possi¬ 
ble trend towards higher Ha EWs at higher redshifts, which 
is extremely important for estimating the sSFR at 2 > 5 
(see also Labbe et al.|2013 de Barros et al.|2014 Smit et al. 


2014[ Smit et al.||2015 l. Independent evidence for high-EW 

nebular lines having a large impact on the broadband fluxes 
of 2 > 3 galaxies was obtained from early near-infrared. 


multi-object spectroscopic campaigns (e.g., Schenker et al. 
20Tl|Holden et al. 2015|). 


The highest-redshift window providing us a largely 
clean measurement of the Ha EW is the redshift range 
2 ~ 5.10-5.40. Here the flux excess due to the redshifted 
nebular emission lines gives rise to significantly redder [3.6] 
— [4.5] colors over this range, making it possible to quantify 
the EW of Ha at 2 > 5.1 in a similar way to that possible in 
the redshift range 2 = 3.8 to 5.0. By examining the [3.6] — 


[4.5] color of galaxies from the Bouwens et al. (20151 catalog 
(see §3.1) over the Great Observatories Origins Deep Survey 
(GOODS)-North and South Fields (Giavalisco et al. 20041 
in this highest redshift window 2 ~ 5.1-5.4 where Ha can be 


cleanly measured^ we can derive approximate constraints 
on the Ha flux and EW at the highest redshift currently 
accessible at reasonable S/N with current facilities. This al¬ 
lows us to obtain better constraints on the evolution of the 
mean Ha EW and sSFR as a function of redshift. 

This paper is structured as follows. In §2, we describe 
the observational data sets and our photometric selection 
of sources in the narrow redshift range. In §3, we examine 
our selected sample of galaxies. We describe the assump¬ 
tions made in deriving the EWs and sSFRs. Finally, in §4, 
we discuss our results and give a summary. We refer to the 
HST F435W, F606W, F775W, F814W, F850LP, F105W, 
F125W, and F160W bands as B435, Veoe, *775, hi 4 , 2350, 
Flos, J 125 and i7i60, respectively. For consistency with pre¬ 
vious work, we adopt the concordance model with Qm ~ 
0.3, Da = 0.7 and Hq = 70 km s“^ Mpc“^. Throughout we 
assume a Salpeter (19551 initial mass function (IMF) be¬ 


tween 0.1—100 Mq. All magnitudes are quoted in the AB 
photometric system (jOke &: Gunn|[T983|). 


2 OBSERVATIONS 
2.1 Data 

In order to select sources in the redshift range 2 ~ 5.1- 
5.4, we make use of the deep optical/ACS and near- 
infrared/WFG3/IR observations over the GOODS-North 
and GOODS-South Fields from three significant Hubble 
Space Telescope (HST) programs: GOODS, ERS (Early Re¬ 
lease Science: Windhorst et al.|2011 |, and CANDELS (Cos¬ 
mic Assembly Near-infrared Deep Extragalactic Legacy Sur¬ 
vey: Grogin et al.|20lT Koekemoer et al.|2011 1. The moder¬ 
ately deep regions over the CANDELS GOODS-North and 
South reach a 5a depth of ~ 27.5 mag in the Fios, J 125 and 
H 160 filters with the HST and cover ~ 100 arcmin^. The deep 
regions over the CANDELS GOODS-North and South reach 
a 5cr depth of ~ 28.5 in the Fios, J 125 an d iLieo bands, cov¬ 
ering ~ 125 arcmin^ (Grogin et al. |201l[ ). HST observations 
with the ACS (Advanced Camera for Surveys) are available 
in the B 435 , V eoe, irrs, hn and 235 0 bands, up to ~ 29 mag 
at 5a in 7 gi 4 ( Bouwens et al.|2015 l. Over the northern ~40 
arcmin^ section of GOODS South ( Windhorst et aL||2011 |, 
deep near-IR observations are available (~28 mag at 5a) in 
the logs, J 125 , and Tlioo bands and also in the R 435 , Vgoe, 
7775 and 2850 bands with ACS. The observations are PSF 
matched to the Hieo band before measuring the colors in 
scalable [Kron |( p80| apertures. 

Essential to the analysis we perform is the 3.6 pm 
and 4.5 pm IRAG observations from the Spitzer Space 
Telescope. Spitzer/TRAC data is from the original GOODS 
program, the Spitzer Extented Deep Survey (SEDS: [Ashbyj 
et al. 2013|), the Spitzer Very Deep Survey (S-CANDELS: 
Ashby et al.|201^ Exploratio n Science Project, the lUDFlO 
program ( [Labbe et al. 2015p , and other programs (such as 
PID10076, PI: Oesch). The Spitzer/IRAC reductions we uti- 


^ At 2 > 5.5, constraints on the Ha EWs are also possible from 
Spitzer/IRAC observations, but would need to rely on the stack¬ 
ing the fluxes of 2 > 5.5 galaxies in the much less sensitive 5.8 
pm band. 
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Figure 1. Left panel: Model spectral energy distribution (blue) and the HST flux measurements (black dots) for one of the sources in 
our selection fitted with the photometric redshift code EAZY ( [Brammer et al.|2008^ . Displayed in the upper left corner: the cumulative 
probability distribution P(z) to find the galaxies at a certain redshift. We only use candidates with P(5.1 < z < 5.4) > 0.85. The 
Spitzer/IRAC bands are not included in the fitting to prevent bias in the [3.6] — [4.5] colors. Right panel: Photometric redshift as 
determined by EAZY against spectroscopic redshift for the sources in [Vanzella et aT)] | |2009] l and [Stark et al.] | |2013[ |. At z < 5, EAZY 
seems to overestimate the redshift by Az/(1 + z) ~ 0.02, while the redshift is underestimated by Az/(1 + z:) > 0.03 for z > 6.0. Though 
we need a far larger sample to confirm the presence of such an offset, the redshift window that we are considering is seemingly the least 
affected. 


lize were generated by [Labbe et al. (2015)| and feature a 
1.8"-diameter FWHM for the PSF. 

Deblending neighboring galaxies in the IRAC ob¬ 
servations and PSF correctons are performed using the 
MOPHONGO software (Labbe et al. 2010a 2010b 2013 20151. 
HST F160W images are used as a high resolution prior to 
construct a model for the contaminating sources, while leav¬ 
ing the normalization of the sources as a free parameter. The 
fluxes of all sources in a radius of 13” are then simultane¬ 
ously fit to best match the IRAC image. Photometry is then 
performed within a 2.0” diameter aperture. The deep IRAC 
imaging from S-CANDELS reaches ~26.8 mag at 5cr in the 
3.6 pm band. 

For 2 = 5.1-5.4 sources not included in the photometric 
catalogs of Bouwens et al. (2015: due to their being located 
in areas of the GOODS fields without B 435 or Fogs/Lios- 
band observations), we made use of the HST photometry 
from the 3D-HST GOODS-North or GOODS-South cata¬ 
logs ( Skelton et al.|[^014 |. This is relevant for 5 2 = 5.1-5.4 
galaxies in our final sample. We refer to Skelton et al. (20141 
for a detailed description. 


2.2 Photometric Redshift Selection 

The Lyman-break selection criteria applied for sources at 2 
~ 5 are as follows: 

(Lfios ~ *775 > 1-2) A (2850 ~ W'leo < 1-3) A 
(Veoe — *775 > 0.8(2850 — Rieo) + 1-2) 

where A denotes the AND symbol. For non-detections, the 
la upper limit is taken as the flux in the dropout band. The 
aforementioned criteria enable to select sources in the range 
2 ~ 4.5 — 5.5. Therefore, a high-redshift boundary is set by 


excluding sources which satisfy the selection criteria for 2 
~ 6 selection (e.g., Bouwens et al. 2015a). Contamination 
from sources at lower redshifts is reduced by requiring that 
the 2 ~ 5 sources have a non-detection (< 2a) in the 5435 
band. Furthermore, we exclude point sources by requiring 
the SExtractor stellarity index to be less that 0.9 (where 
0 and 1 correspond to extended and point sources, respec¬ 
tively). Utilizing these selection criteria results in an initial 
U-drop sample of 1567 sources ( Bouwens et al.|[2015 K 

In order to be able to measure the Ha EW for 2 > 5 
galaxies, it is easiest to only make use of galaxies in the 
narrow redshift window 2 = 5.10-5.40. We therefore use 
photometric redshifts to identify a sub-sample of galaxies 
in this window. The photometric redshifts for our sample 
are determined using the EAZY photometric redshift code 
(Brammer et al. 20081, which compares photometric data 
with synthetic photometry of galaxies for various template 
spectra and redshift ranges. The best-fit redshift is then de¬ 
rived from a statistical analysis of the differences between 
both data sets. The aforementioned 8 HST bands (B435, 
Veoe, 7775, 7814, 2850 , Yio 5 , 7i25, Rioo) are used to derive the 
best-fit photmetric redshifts. The IRAC photometry is ex¬ 
cluded from this fitting to avoid introducing any bias in the 
measured [3.6] — [4.5] color. This reduces the sample size to 
393 sources. 


Figurej^shows an example output of EAZY for a source 
at 2 ~ 5.2 and a comparison of our estimated EAZY photo¬ 
metric redshifts for a sample of 2 = 4 — 6 spectroscopically 
confirmed sources from Stark et al. (2010: see D. Stark et 
al. 2015, in prep) and Vanzella et al. (2009). The scatter 
around the one-to-one relation is A 2 /(l-|- 2 ) = 0.036, which 
provides confidence that we can select sources in the narrow 
redshift range 2 = 5.1-5.4. 
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Figure 2. Nebular emission line contamination in the 3.6 /im 
and 4.5 /i.m photometric filters. Top panel : The redshift ranges 
over which the dominant nebular emission lines, Ho;, H/3, [OIII], 
[Nil], and [SII], contribute to the 3.6 pm and 4.5 pm flux mea¬ 
surements. Lower panel : The predicted [3.6] — [4.5] color due 
to various nebular emission lines as a function of redshift. The 
solid black circles indicate sources that are selected in the red¬ 
shift range Zp^ot = 5.1-5.4, where Ha lies in the 4.5 pm filter, 
while 3.6 pm is devoid of strong nebular emission. Their observed 
[3.6] — [4.5] colors are primarily very red. The open circles are 
the colors for the spectroscopic sample found in [Vanzella et al.| 
l |2009[ |, which we use to estimate a stellar continuum color of 
~ 0.00 mag. The dotted line indicates the expected color for a 
strong evolution of the r est-frame equivalent wi dth, EWo(Ho) oc 
(l+z)^ ® A, according to Fumagalli et al. 1 2012 I. Our selection of 
sources in the redshift range z = 5.1-5.4 has a mean [3.6] — [4.5] 
color of 0.31 it 0.05 mag and 0.35 it 0.07 mag for our photometric- 
redshift and spectroscopic samples, respectively, implying a mean 


EW(Ha-|-[NII]-|-[SII]) of 665 it 53 A and 707 it 74 A, respectively, 
for sources in these samples (and 638 it 118 A and 855 it 179 A, 
respectively, using direct SED fits). Four selected sources show 
even far redder colors than predicted with the strong evolution 
model. The four reddest sources have a mean color of 0.66 it 0.06 
mag, with a notably high EW of 1743 it 221 A. 


We select very bright sources with the requirements that 
S/N(J7 i6o) > 7 a S(7yi6o)/N(3.6 /rm) > 3 A S(7/i6o)/N(4.5 
/tm) > 3, where S/N is the signal to noise ratio. Our selection 
of sources based on their measured flux in the Tfieo band 
and measured noise in the Spitzer/IRAC bands allows us 
to include sources in our analysis which we would expect to 
show up prominently in the Spitzer/IRAC bands. Basing the 
selection on the measured flux in the Spitzer/IRAC bands 
would bias our measurement of the colors. We also discard 
41 sources for which the contamination by nearby objects is 



Figure 3. [3.6] — [4.5] color distribution for sources in our z = 5.1- 
5.4 primary selection (red histogram) and z = 4.4-5.0 control sam¬ 
ple (blue histogram: Appendix A). Sources where Hq!-|-[NII] is in 
the [3.6] band uniformly have moderately blue colors, whereas 
sources with Hq:-|-[NII] in the [4.5] band uniformly have moder¬ 
ately red colors. The impact of the Hq:-|-[NII] emission lines on 
the IRAC colors is quite clear. Deriving the Hq!-|-[NII] + [SII] EW 
by contrasting the observed colors of sources in the two samples 
should produce a very robust result. 


higher than 200 njy or sources with a poor IRAC deblending 
(X^ parameter less than 0.2), reducing the sample to 101 
bright sources. 

Finally, we require that the EAZY redshift probability 
distribution prefers a redshift in the range z = 5.1-5.4 at 
> 85%, providing a sample with 11 sources in the redshift 
range z ~ 5.1-5.4 presented in Table 


2.3 Spectroscopic Redshift Selection 


In addition to making use of sources very likely to lie in 
the redshift range z = 5.1-5.4 using our photometry for the 
sources and the redshift likelihood distributions we derive, 
we can also make use of sources known to lie in the redshift 
range z = 5.1-5.4 from available spectroscopy (D. Stark et 
al. 2015, in prep). Using spectroscopic redshifts, we can be 
even more certain that the sources we are using lie in the 
narrow redshift range z = 5.10-5.40 required for our desired 
measurement of the Ha flux. 

One potential drawback to the inclusion of such sources 
in the present study is that we might be working with a bi¬ 
ased sample, given that essentially all of the spectroscopic 
redshift measurements we utilize come from Lya, and it is 
not clear a priori that the study of such a sample might 
bias the mean Ha EW we measure to higher values. Fortu¬ 
nately, as we show in a separate study (Smit et al. 2015, in 
preparation), the mean Ha-|-[NH]-|-[SH] EW measured for 
both photometric-redshift and spectroscopic samples are es¬ 
sentially identical. 

Cross-correlating the source catalogs of|Bouwens et ah] 


(20151 and Skelton et al. (20141 with the spectroscopic cat- 
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alog of D. Stark et al. (2015, in prep) and Vanzella et al. 
(20091, we identified 13 2 ; = 5.10 — 5.40 galaxies that we 


can use for our study, with one source overlapping with our 
photometric selection (see Table [^. As in our photometric 
redshift selection (§2.2), we exclude sources where flux from 
neighboring sources significantly contaminate the photomet¬ 
ric apertures for our 2 = 5.10-5.40 sample (3 sources). 10 of 
the sources in the desired redshift range were from D. Stark 
et al. (2015, in prep) redshift compilation, while 3 came from 
the Vanzella et al. (2009) compilation. 

We will make use of sources from both our high-quality 
photometric redshift sample and spectroscopic sample for 
the analyses that follow. 


3 RESULTS 

Using the selection discussed in §2.2 and §2.3 we have iso¬ 
lated a sample of 21 galaxies that have redshifts in the nar¬ 
row redshift range 2 ~ 5.1 — 5.4 at high conhdence. In this 
section we will discuss the IRAC [3.6] — [4.5] colors of these 
galaxies and compare them with the IRAC colors of a spec¬ 
troscopically conhrmed sample of 2 ~ 4.5 galaxies. 


3.1 Mean [3.6] — [4.5] Color for 2 = 5.1 — 5.4 
Galaxies 

Our selection of galaxies in the redshift range 2 ~ 5.1 — 
5.4 allows us to solve for the rest-frame EW of the nebular 
emission lines using the [3.6] — ]4.5] color. Assuming Fv « 
constant (i.e. Fx oc A“^) we approximate the observed flux 
Fv,ohB, in the IRAC filters, for a given EW, by 

Fi/ ohs — Ey continuum ' iT 


Sew ~ ( 1 

\ lines,i 


EWo,i • (1 -l- z) • R(Aobs,i) 
Aobs.t f R,{X)/\dX 


( 1 ) 


where Eh,continuum is the stellar continuum flux and 
Aobs.i is the observed wavelength of the nebular emission 
lines (Smit et al 1014| ). R{\) denotes the response curve of 
the Alter. We use all nebular emission lines tabulated in 
Anders & Fritze-v. Alvensleben ( 2003| ) and the hydrogen 
Balmer lines for the modeling of the [3.6]— ]4.5] color. We 
fix the line intensities relative to H/3 according to the val¬ 
ues in [Anders fc Fritze-v. Alvensleben] ( |2003| ) for sub-solar 
metallicity 0.2 Zq and assuming case B recombination. The 
observed ]3.6] — ]4.5] color can then be modeled by 


[3.6]-[4.5] = ([3.6]-[4.5])co„ti: 


.-2.51ogio (^) . (2) 

V®4.5 / 

Figure [^illustrates the contamination by nebular emis¬ 
sion lines in the photometric filters as a function of redshift. 
The top panel shows the redshift ranges where the dominant 
emission lines, Ha, H/3, [OHI]AA4959, 5007, [Nil], and [SH] 
fall in the IRAC 3.6 /rm and 4.5 /rm passbands. The lower 
panel shows the expected [3.6] — [4.5] color in the presence 
of nebular line emission. The dotted line illustrates the case 
when the rest frame EW evolves strongly with redshift. Con- 


oc ( 1 - 1 - 2 )^'® A for star forming galaxies at 0 < 2 < 2. Galax¬ 
ies are expected to become quite red in the redshift range 2 
= 5.1 — 5.4, due to contamination of the 4.5 /rm flux by the 
Ha line and no strong nebular emission line contamination 
in the 3.6 fim band. 

The selected galaxies in the redshift range 2 ~ 5.1 — 5.4 
allow us to deduce the EW of the nebular emission lines in 2 
> 5 objects from the fluxes in the 3.6 fj,m and 4.5 fj,m bands. 
We observe a mean [3.6] — [4.5] color of 0.31 ± 0.05 mag 
for the selected sources in our photometric-redshift sample 
and a mean [3.6] — [4.5] color of 0.35 ± 0.07 mag for sources 
in our spectroscopic sample. In both cases, we estimate the 
uncertainty here through bootstrap resampling. 

We estimate the [3.6] — [4.5] continuum color by con¬ 
trasting the mean [3.6] — [4.5] colors observed for 2 = 4.4- 
5.0 spectroscopically-confirmed sources from [Vanzella et al.j 


(20091, Shim et al. (20111, and 


Stark et al. 


(20131, i.e.. 


—0.32 ±0.03 mag, with the mean [3.6] —[4.5] colors observed 
for our selected 2 ~ 5.1 — 5.4 sample and attribute any dif¬ 
ferences in the [3.6]— [4.5] colors to the impact of flux from 
the Ha ± [Nil] ± [SH] lines. We estimate a color for the 
stellar continuum of ]3.6] — [4.5] of 0.00 ± 0.04 mag. See 
Figure [^ for an illustration of the differences between the 
two samples. 

Four sources in our selection show considerably redder 
[3.6] — [4.5] colors than we can infer from our model SED 
with strong line emission, i.e., 0.66 ± 0.06 mag. The very 
red [3.6] — [4.5] colors observed for many sources in our 
2 = 5.1 — 5.4 sample suggests that the IRAC [3.6] — [4.5] 
color may provide significant leverage in terms of identifying 
galaxies that specifically lie in the narrow redshift range 2 
~ 5.1 — 5.4. This method has previously been explored at 
2 ~ 7 — 8 by [Smit et al.| ( [2015 1 and [Roberts-Borsani et al.| 
(20151 based on the impact of [OHI] on the IRAC fluxes at 


these high redshifts. 


3.2 Mean Ha±[NII] + [SII] EW of 2 ~ 5 Galaxies 

By comparing the mean [3.6] — ]4.5] color of our 2 = 5.10- 
5.40 sample with the mean [3.6] — [4.5] color of our 2 = 4.4- 
5.0 sample, we find an overall color difference of 0.68T0.08 
mag relative to our 2 = 5.1-5.4 photometric sample and 
0.69T0.09 mag relative to our 2 = 5.1-5.4 spectroscopic sam¬ 
ple. Comparing the color difference we observe with that 
predicted based on simple model s pectra with an Ha EW 
of 595 A and line ratios set by the Anders & Fritze-v. Al 


vensleben (20031 model {dotted line in Figure 2), we infer 


sistent with Fumagalli et al. (2012) who And that EWo(Ha 


an approximate Ha EW of 557 ± 44 A for our 2 = 5.1-5.4 
photometric sample and 592 ± 62 A for our 2 = 5.1-5.4 
spectroscopic sample. 

Given that flux in the [Nil] and [SH] lines would add 
to the color difference observed between our 2 = 4.4-5.0 
and 2 = 5.1-5.4 samples and cannot be determined sepa¬ 
rately, it is best to quote a constraint on the mean EW of 
Ha± [Nil] ±[SH]. I Anders fc Fritze-v. Alvensleben] ( [2003[ ) pre¬ 
dict a contribution of 6.8% from [Nil] and 9.5% from ]SH]. 
This is in good agreement with an observed ratio of [Nil] /Ha 
of 0.05-0.09 in 2 ~ 2.3 galaxies with stellar masses in the 
range log{Mt/M q) = 9.15 — 9.94 by [Sanders et al.1[2015| . If 
we correct for this, the mean EW in Ha±[NH]±[SH] is 665 
± 53 A for our photometric sample and 707 ± 74 A for our 
spectroscopic sample. 
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Table 1. Our selections of sources with a high probability of lying in the redshift range 5.1 < 2 : < 5.4 from the photometry.^ 








[3.6] 

— 

[4.5] 

Ha+[NII] + [SII] 

sSFR 

logicM 




ID 

RA 

DEC 

^phot 

(ma^ 

;) 

EW [A]'^ 

[Gyr-i] 

[Mq] 

Kv 


GNDV-7133823953 

12:37:13.38 

62:12:39.5 

5.2 

± 

0.1 

0.5 

± 

0.3 

672 ± 505 

30tl° 

8.84l0;ll 

-21.3 

± 

0.1 

GNDV-7128013231 

12:37:12.80 

62:11:32.3 

5.3 

± 

0.2 

0.1 

± 

0.2 

190 ± 254 


9.24l0;l° 

-20.8 

± 

0.1 

GNDV-7033233179‘^ 

12:37:03.32 

62:13:31.8 

5.3 

± 

0.1 

0.4 

± 

0.2 

719 ± 368 

24l?i 

8.93l°;l° 

-20.8 

± 

0.1 

GNDV-6302234526 

12:36:30.22 

62:13:45.3 

5.2 

± 

0.1 

0.3 

± 

0.1 

385 ±170 

n+6 

9.38l°i^ 

-20.6 

± 

0.1 

GNDV-6285841077 

12:36:28.58 

62:14:10.8 

5.3 

± 

0.1 

0.2 

± 

0.2 

557 ± 449 


Q 7’c*-|-0.32 
'’^-0.20 

-20.7 

± 

0.1 

GNWV-6514085687 

12:36:51.40 

62:08:56.9 

5.3 

± 

0.1 

0.0 

± 

0.1 

25 ± 167 

111 

iU.lZ_Q QQ 

-20.6 

± 

0.1 

GNWV-6121502518 

12:36:12.15 

62:10:25.2 

5.3 

± 

0.1 

0.4 

± 

0.3 

917 ± 573 

^+5 

'^-10 

Q 2Q+0-25 

y.zy_o 62 

-20.8 

± 

0.1 

GNWV-6095211615 

12:36:09.52 

62:11:16.2 

5.2 

± 

0.1 

0.4 

± 

0.1 

913 ± 274 

lQ-hl2 

Q 24+0-27 
^•^^-0.20 

-21.6 

± 

0.1 

GNDV-3756634257 

12:37:05.66 

62:13:42.6 

5.3 

± 

0.1 

0.1 

± 

0.1 

274 ± 180 

^D_19 

9.55l°i| 

-21.2 

± 

0.1 

GNDV-6325033158 

12:36:32.50 

62:13:31.6 

5.3 

± 

0.1 

0.4 

± 

0.1 

882 ± 262 

9l| 

8.79l0;§^ 

-20.9 

± 

0.1 

GSDV-2332672480 

03:32:33.26 

-27:47:24.8 

5.3 

± 

0.1 

0.6 

± 

0.1 

1480 ± 110 

111 

10.06l°;°« 

-20.6 

± 

0.1 


^ To identify those sources with the highest probability of lying in the redshift range 2 = 5.10 — 5.40, sources 
are required to have P(5.1 < 2 < 5.4) > 0.85. 

^ The 2850 band magnitude is used to derive the intrinsic UV luminosity. 

^ The estimated EW for individual sources is derived by comparing the 4.5 |i.m flux with that derived from 
FAST (excluding the 4.5 |i.m flux from the fits). 

Spectroscopically confirmed to be at 2 = 5.21 (D. Stark et al. 2015, in prep). 


Table 2. Sample of spectroscopically confirmed sources in the redshift range 5.1 < 2 < 5.4.^ 


ID 

RA 

DEC 

'^spec 

^phot 

[3.6] - 
(maj 

[4.5] 

s) 

Ha±[NII[±[SII[ 
EW [A]** 

sSFR 

[Gyr-i] 

logic Af 

[Mq] 

Muv 


GNDV-6554953313 

12:36:55.49 

62:15:33.1 

5.19 

4.9 

± 

0.1 

0.6 

± 

0.2 

1579 ± 359 

lot® 

y-'jy_0.28 

-21.3 

± 

0.1 

GNDV-7033233179 

12:37:03.32 

62:13:31.8 

5.21 

5.3 

± 

0.1 

0.4 

± 

0.2 

719 ± 368 

^^-11 

0 QO+0.10 

^•^'^-0.17 

-20.8 

± 

0.1 

GNDV-7027322916 

12:37:02.73 

62:12:29.2 

5.23 

5.5 

± 

0.1 

0.5 

± 

0.1 

1000 ± 257 

07+12 


-20.6 

± 

0.1 

GNDV-6375223629 

12:36:37.52 

62:12:36.3 

5.18 

5.4 

± 

0.1 

0.0 

± 

0.3 

<419== 

21 + 19 

8.62^07^ 

-20.0 

± 

0.1 

GNDV-6553954912 

12:36:55.39 

62:15:49.1 

5.19 

5.1 

± 

0.2 

0.8 

± 

0.4 

2568 ± 1561 

7+9 

'-5 

8.29tg;g 

-19.8 

± 

0.1 

GNWV-7347782930 

12:37:34.77 

62:18:29.3 

5.32 

5.2 

± 

0.2 

0.6 

± 

0.3 

1344 ± 666 


8.90tgl^ 

-19.9 

± 

0.1 

ERSV-2213040511 

03:32:21.30 

-27:40:51.2 

5.29 

5.3 

± 

0.1 

0.5 

± 

0.2 

1075 ± 394 

l+lo 

9 11+0.27 

-20.9 

± 

0.1 

GSWV-2454254386 

03:32:45.43 

-27:54:38.6 

5.38 

5.4 

± 

0.1 

0.2 

± 

0.1 

289 ± 75 


8-96to.66 

-21.6 

± 

0.1 

GND6418‘^ 

12:36:18.19 

62:10:21.9 

5.28 

5.6 

± 

0.1 

0.3 

± 

0.1 

309 ± 88 


9 99+0-20 

y-yy_o.oi 

-21.6 

± 

0.1 

GND33928 

12:37:36.87 

62:18:55.9 

5.35 

5.2 

± 

0.1 

0.1 

± 

0.2 

<222== 

19+17 

1^-16 

9 09+®'^^ 

-20.4 

± 

0.1 

GND29175 

12:37:31.45 

62:17:08.3 

5.25 

5.3 

± 

0.2 

0.3 

± 

0.4 

666 ± 929 

3+7 


-19.9 

± 

0.1 

GND12038 

12:36:26.49 

62:12:07.4 

5.20 

5.2 

± 

0.4 

0.0 

± 

0.2 

254 ± 384 

lOt^" 

9 11+0-43 

-20.3 

± 

0.1 

GS48361'i 

03:32:16.55 

-27:41:03.2 

5.25 

5.5 

± 

0.2 

0.3 

± 

0.2 

668 ± 370 

6t® 

y.oz,_o 53 

-20.7 

± 

0.1 


^ The spectroscopic redshifts are obtained by cross-correlating the |Bouwens et al.| ( |20f^ and [Skelton et aT] 
( |2014[ | catalogs with the spectroscopic catalog of D. Stark et al. (2015, in prep) and Vanzella et al. (2009). 

^ The estimated EW for individual sources is derived by comparing the 4.5 pm flux with that derived from 
FAST (excluding the 4.5 pm flux from the fits). 

^ We give the error in the Hq;+[NII] + [SII] EW as a Icr upper limit when the inferred value is negative. 

Spectroscopic redshift measurement is based on the identification of a probable absorption line and hence 
less confident than the other spectroscopic redshift measurements included in this table (D. Stark et al. 2015, 
in prep; Vanzella et al. 2009). 


We can also derive Ha+[NII] + [SII] EWs for individual 
sources in our photometric and spectroscopic 2 : = 5.1-5.4 
samples. In computing the EWs for individual sources, we 
use FAST to fit the observed SEDs of individual sources ex¬ 
cluding the 4.5 pm band which is contaminated by Ha emis¬ 
sion. Then, by comparing the observed 4.5 pm flux with the 
expected 4.5 pm flux (without including emission lines in 


the FAST modeling), we derive EWs for individual sources. 
The results are presented in Tables 1 and 2. The mean 
HQ;-|-[NII]-t-[SII] EW we derive for our photometric sample 
is 638 ± 118 A, while we find 855 ± 179 A for our spectro¬ 
scopic sample. If we follow the model results from jAndersj 
& Fritze-v. Alvensleben ( 2003| > and suppose that 16.3% of 


the 4.5 pm excess derives from [Nil] and [SII], the excesses 


© 2015 RAS, MNRAS 000,[iyiTl 





























Ha EW at z ^ 5 from IRAC Colors 7 



1000 


100 


10 


1111111 

1 1 1 1 1 1 1 

" A 

^ ' 

9-^ 

] 

i 






: & 

■ 

- / 

• This work 

/ 

□ Stark et al. 2013 

- / 

<> Schenker et al. 2013 

/ 

A Shim et al. 2011 


X Erb et al. 2006 


O Fumaqalli et al. 2012 

_1_1_1_1_1_1_1_L_ 

-Fumagalli et al. 2012 


0 2 4 6 8 

Redshift 


Figure 4. Hq;+[NII]+[SII] EWs as function of redshift derived from the mean [3.6] — [4.5] color of our source selection. Several estimates 
from the literature are indicated for reference < |Erb e t al.||2006| |Shim et al.||2011| |Fumagalli et al.||2012| |Schenker et al.||2013|Stark et| 
|al.|2013| |. The upper and lower EW determination for Stark et al. (2013) excludes and does not exclude the contaminated IRAC band 
in deriving the stellar continuum required to derive the EW for the Hq:+[NII]+[SII] line. We assume line ratios as listed by [Anders &| 
[Fritze -V. Alvensleben| l |2003[ ) and that the stellar continuum of our z = 5.10-5.40 sample has a [3.6] — [4.5] color of 0.00 ±0.04 magnitude. 
Redder [3.6] — [4.5] colors will therefore be due to the emission lines contaminating the 4.5 fim flux. The measured Ha±[NII]±[SII] EW 
shown here represents the weighted mean of our EW estimates from our photometric and spectroscopic z = 5.1-5.4 samples (i.e., 684±51 
A) and is higher than values derived at lower redshifts, suggesting stronger line emission at z ~ 5. The evolution of Ha±[NII]±[SII] EW 
for the indicated stellar mass range found by |Fumagalli et al.| (|2012|| is extrapolated (given by the dashed line) and is consistent with 
our inferred H«±[NII]±[SII] EW. 


we derive suggest Ha EWs of 534 ± 99 A and 715 ± 150 A, 
respectively. 

The four reddest sources have a mean [3.6] — 
[4.5] color of 0.66 ± 0.06 mag, consistent with an 
EWo(Ha±]NII]±[SII]) of 1743 ± 221 A, equivalent to 
EWo(Ha) = 1458 ± 185 A for the above line ratios. 

Figure]^ shows several values for the Ha ± [Nil] ± [SII] 
EWs with redshift from the literature. The black line gives 
the evolutio n of the Ha ± [Nil] ± [SII] derived by FugamalII| 
et al. 120121 for galaxies with masses M ~ 10^° — 10^^“ Mq, 


which we extrapolate to higher redshifts and lower masses. 
Keeping in mind that the EW(Ha ± ]NII] ± [SII]) as func¬ 
tion of the redshift is higher for sources with lower steller 
masses our result is consistent with the extrapolation from 


Fumagalli et al. (20121 and the high EWs derived by Shim 


et al. (20111, Stark et al. (20131 and Schenker et al. (20131. 


Uncertainties in the photometric redshifts for our z = 
5.10-5.40 sample can lead to a systematic underestimate of 
the Ha±[NH]±[SH] flux, if it causes us to include sources 
which lie outside the desired range. For z < 5.1 sources, 
the 3.6 pm band will be contaminated by Ha ± [Nil] ± 
[SII] emission. Meanwhile, for z > 5.4 sources, H/3 ± [OIII] 
emission will contribute to the 3.6pm band. In both cases, 
the [3.6] — [4.5] color will be much bluer, causing us to infer a 
substantially lower EW for that source, than is truly present. 


3.3 Specific Star Formation Rates 


A clean measurement of the stellar continuum emission is 
essential for deriving the sSFR. In our target redshift range 
the flux in the 3.6 pm band can be used for this purpose, 
while the 4.5 pm band is contaminated by emission lines and 
should be left out in stellar population modeling. We derive 
the mean sSFR for our selected sources using similar method 


as described in Stark et al. (20131. We derive stellar masses 
using the modeling code FAST ( Kriek et al. [2009 1, which fits 
stellar population synthesis templates from |Bruzual fc Char-| 
lot (20031 to broadband photometry. For consistency, we as- 


Salpeter (19551 IMF with 0.1—100 Mq, a sub-solar 


metallicity Z = 0.2 Zq and dust attenuation from |Calzetti| 
et al. (20001. The ages range from 10 Myr to the age of 


the universe at the redshift of the source. The star forma¬ 
tion history is assumed to be constant and the dust content 
is varied between Ay = 0 — 3 mag. In order to avoid the 
degeneracy between the ages of the galaxies and their dust 
content in the FAST modeling of the sSFR, we derive star 
formation rates directly from the UV continuum, using the 


Kennicutt (19981 relation and fixing the dust extinction us¬ 
ing the relationship from Meurer et al. (19991 between Ay 


and the UV-continuum slope /3. From this method, we derive 
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log [Moss/Mq] 

Figure 5. The estimated Hq!+[NII] + [SII] equivalent width from the 4.5 fim excess vs. the estimated stellar mass (red squares with Itr 
error bars). The mean Ha+[NII] + [SII] EW we estimate for sources with an estimated stellar mass > 10® ®Mq is 522 it 279 A (blue 
diamond), while the mean Ha+[NII]+[SII] EW we estimate for sources with an estimated stellar mass < is 773 it 136 A (green 

triangle). The present results provide no strong evidence for a correlation between the inferred Ha+[NII]+[SII] EWs and the stellar mass. 
However, we emphasize that this may change as the samples and dynamical masses become larger. 


a median sSFR of I 7 I 5 Glyr ^ (individual sSFR estimates 
are listed in Table and Table |^. 

We also explore the use of a SMC dust-law, which might 
seem more appropriate for Lyman-break galaxies at z ~ 5 —6 
based on recent ALMA results by|Capak et al.| (|2015[ ). We 
fit FAST models as above with a Noll et al. (20091 dust law 
with parameters Eb = 0.01 and S = —0.42 and find a sSFR 
of ~ 13 Gyr“^. 

It is useful to compare the median sSFR estimates we 
find here, 17lg Gyr“^ and 13 Gyr“^, with previous esti¬ 
mates. This represents a 7^.2 ^ increa se in sSFR compared 
to the 2.4 Gyr”'^ value found at z ~ 2 ( Reddy et al.|2012a l, 
supporting a significant evolution in the sSFR. The typical 
stellar mass for galaxies in our 2 : = 5.1-5.4 selection is ~10® 
Mq , so we will make our comparison with previous measures 
at this stellar mass. Both Stark et al. (2013) and Gonzalez et 
al. (2014) find that the typical galaxy with this stellar mass 
has a UV luminosity Muv of —20 mag. Accounting for a 
factor of ~2 mean dust attenuation at this luminosity (as 
implied by the Bouwens et al. 2014 ~ —1.9 results), the 

equivalent sSFR is ~11 Gyr“^. Somewhat similarly, Salmon 
et al. (2015) derive a sSFR of Gyr“^, which is again 

somewhat lower than we find here. 


All things being equal, we would expect the sSFR es¬ 
timates we derive here to be more accurate than previous 
estimates, given our precise knowledge of the redshifts and 
hence the position of nebular emission lines within galaxies. 
The impact of the lines on the mass and sSFR estimates 
could be as large as a factor of ~ 1.5, allowing for an approx¬ 
imate reconciliation of the present sSFR estimates with that 
from previous work. However, the present sample of 2 ~ 5.1- 


5.4 galaxies is still quite small, and therefore expansion of 
the present sample would certainly be helpful for improving 
our sSFR estimate. 


3.4 Possible Dependence of the EW of 
Ha+[NII]-|-[SII] on the Stellar Mass 

By fitting to the photometry of all passbands uncontami¬ 
nated by the strong Ha-|-[NII]-|-[SII] nebular emission lines, 
we can estimate stellar masses for sources in our 2 = 5.1- 

5.4 samples, as described in the previous section. As these 
sources are distributed over a wide range of stellar mass, i.e., 
10 ®-® Mq to 10^° '* Mq, we can go beyond a simple determi¬ 
nation of the mean Ha-|-[NII]-|-[SII] EW for 2 ~ 5 galaxies 
and look at whether there is a dependence on stellar mass. 

Any significant dependence on stellar mass would be 
noteworthy, as it could point to a significant mass or scale 
dependence to the star formation histories of galaxies. While 
such a scale dependence could be expected if there are sig¬ 
nificant feedback effects at early times (e.g., Bowler et al. 
2014; Bouwens et al. 2015), many simulations (e.g., Finla- 
tor et al. 2011 ) predict that galaxies build up their stellar 
mass in a relatively self-similar manner, independent on the 
overall stellar mass. 

We examine the evidence for such a correlation between 
the Hq-|-[NII]-|-[SII] EW and the stellar mass in Figure]^ 
plotting the Ha-|-[NII]-|-[SII] EWs we estimate for individ¬ 
ual sources against the stellar masses we estimate for the 
same sources. Also shown on this figure is the average Ha 
EW we measure for galaxies with estimated stellar masses 
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<10® ® Mq (773 ± 136 A) and for those with estimated 
stellar masses >10®'® Mq (522 ± 279 A). 

At face value, these results do not provide any strong 
evidence (< 2 ct) for a correlation between the Ha EW and 
the stellar mass. While it is possible that a slight correlation 
might be expected (i.e., since massive galaxies would be the 
first to experience a slowing in their growth rate due to 
feedback-type effects), no strong trend is evident. A similar 
conclusion can be drawn by fitting Ha-|-[NII]-|-[SII] EWs vs. 
stellar mass relation to a straight line. 


4 DISCUSSION AND SUMMARY 

In this paper, we derive the Ha-|-[NII]-|-[SII] EW and the 
sSFR by selecting galaxies at « ~ 5.1-5.4. In doing so, we 
make use of the highest redshift window allowing for a clean 
measurement of both the Ha flux and the stellar continuum. 
Reliable estimates of the stellar mass and sSFR require a 
clean measurement of the stellar continuum. At z — 5.4- 
6.6, rest-frame optical flux information on galaxies - as de¬ 
rived from the Spitzer/IRAC data - can be quite ambiguous 
to interpret, due to a significant contribution from nebular 
emission lines (Ha, H(?, [OIII], [Nil], [SH]) in both sensitive 
Spitzer/IRAC channels (i.e., [3.6] and [4.5]), and at 2 > 6.6, 
the Ha line redshifts out of the [4.5] channel. 

Observing galaxies in the 2 = 5.1-5.4 redshift interval is 
useful as a contrast and for interpreting the observations of 
galaxies in the redshift interval 2 = 3.8-5.0. In both intervals, 
the Spitzer/IRAC fluxes are expected to be dominated by 
a stellar continuum contribution and a contribution from 
Ha. However, the Ha line will contribute to the measured 
flux in a different Spitzer/IRAC band at 2 = 5.1-5.4 ([4.5]) 
than at 2 = 3.8-5.0 ([3.6]), so the colors of galaxies in the two 
samples can be contrasted and used to set strong constraints 
on the rest-frame EW of Ha. 

In our utilization of the 2 = 5.1-5.4 redshift interval 
to study the Ha EW and sSFR, we have selected 11 bright 
sources over the CANDLES GOODS-North and GOODS- 
South Fields satisfying the criteria S(i7i6o)/N(3.6 /rm) > 
3, S(Hi6o)/N(4.5 /rm) > 3, (Veoe — *775 > 1.2 ), (2850 — 
H 160 < 1.3), (Veoe — *775 > 0.8 (2350 — Hieo) + 1-2), and 
S/N(R 435 ) < 2. The candidates are required to have > 85% 
probability of lying in the redshift range 2 ~ 5.1 — 5.4. We 
have supplemented this sample with 13 2 = 5.10 — 5.40 
sources from the spectroscopic redshift sample of Vanzella 
et al. (2009) and D. Stark et al. (2015, in prep). 

We find a mean rest-frame EW(Ha-|-[NH[-|-[SH]) of 665 
± 53 A for our photometric sample and 707 ± 74 A for our 
spectroscopic sample based on the mean [3.6] — [4.5] colors 
of these samples. Assuming that 84% of the Ha-|-[NH]-|-[SH[ 
line flux is Ha, we further derive a Ha EW of ~ 557±44 
A and 592 ±62 A for our photometric and spectroscopic 
sample, respectively. Our estimate is consistent with the 
( 1 ± 2 )^'® power law derived for a strong line-emitter model 


to the 2.4 Gyr“^ value found at 2 ~ 2 ( [Reddy et al.|2012a l, 
supporting a significant evolution in the sSFR. Our estimate 
is in agreement with the theoretical model of [Neistein fc| 
Dekel (20081, matching the increasing specific inflow rate of 


baryonic particles. 

We emphasize the Ha±[NH]±[SH[ EWs and sSFR we 
derive here for our photometric sample is effectively a lower 
limit on the true value, as the inclusion of any sources in our 
photometric selection outside the desired redshift range due 
to photometric redshift uncertainties would result in a bluer 
mean [3.6] — [4.5] color and higher redshift optical flux for 
the average source. 

We also take advantage of the stellar population model¬ 
ing we do of individual sources in our sample and the range 
of estimated masses to look at a possible correlation between 
the Ha±[NH]±[SH] EW of sources in our selection and their 
stellar masses. We find no strong evidence (<2 ct) for there 
being a correlation between the EW of Ha±[NH[±[SH] EWs 
and the stellar mass. However, we caution that our sample 
sizes and dynamic range are limited, and so a correlation 
may be evident when examining large sample sizes or a wider 
dynamic range. 

More accurate results would require spectroscopic red- 
shifts for a larger number of sources. Though our target red¬ 
shift range allows for a relatively clean measurement, con¬ 
structing large samples of 2 ~ 5.1 — 5.4 galaxies is chal¬ 
lenging due to the narrow redshift window we are consider¬ 
ing. We expect significant progress in the future as a result 
of future samples with the MUSE spectrograph (Bacon et 
al. 2015). Furthermore, ultra-deep Spitzer/IRAC data over 
the GOODS-N and GOODS-S fields will become available 
through the GOODS Re-ionization Era wide-Area Treasury 
from Spitzer (GREATS, PI: Labbe) program (2014). 
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P. R. M., Ferguson H. C., 2006, ApJ, 651, 24 


APPENDIX A: Z = 4.4 - 5.0 REFERENCE 
SAMPLE 

In order to estimate the mean stellar continuum color 
[3.6] — [4.5] for z ~ 5 (§3.1), we make use of a spectroscopic 
selection of sources at z = 4.4-5.0 from Vanzella et al. (2009), 
Shim et al. (2011), and Stark et al. (2013). We obtain a sam¬ 
ple of 30 2 = 4.4-5.0 sources by cross-correlating the source 


catalogs of Bouwens et al. (20151 and Skelton et al. (20141 


with the spectroscopic catalog of Shim et al. (2011), Stark 


et al. (2013) and Vanzella et al. (20091. As in section 2.3 


we exclude any sources which are reported to have detected 
X-ray counterparts or AGN emission lines by [Shim et al.j 
(20111, Vanzella et al. (20091 and D. Stark et al. (2013). The 


excluded sources including those with the following coor¬ 
dinates (03:32:29.29, -27:56:19.5; 03:32:44.11, -27:54:52.5; 
12:36:42.05, 62:13:31.7; 03:32:33.77, -27:52:23.7). We tabu¬ 
late the measured [3.6] — [4.5] colors and coordinates of the 
sources we utilize in Table iati 
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Table Al. Our reference sample of spectroscopically confirmed sources in the redshift range 2 ~ 4.4 — 5.0. 
The sources are obtained by matching the spectroscopic redshift sample listed in|Vanzella et al.|||2009[| w ith 
the|Bouwens et al.| ( [201^ catalog. Source IDs are as in the |Bouwens et al.| < |2015| |, |Shim et al.| ( |2011| |, or 
[Stark et al.| ( |2013| l catalogs. 


ID 

RA 

DEG 

Zspec 

[3.6] - 

- [4.5] 

GSWV-2426242897 

03:32:42.62 

-27:54:28.97 

4.400 

-0.5 

± 

0.1 

GSDV-2228872758 

03:32:22.88 

-27:47:27.58 

4.440 

1 

p 

CO 

± 

0.1 

GSDV-2229762901 

03:32:22.97 

-27:46:29.01 

4.500 

-0.4 

± 

0.1 

ERSV-2285605575 

03:32:28.56 

-27:40:55.75 

4.597 

-0.2 

± 

0.2 

GSDV-2169812296 

03:32:16.98 

-27:51:22.96 

4.600 

-0.4 

± 

0.1 

GSWV-2475822816 

03:32:47.58 

-27:52:28.16 

4.758 

1 

p 

± 

0.2 

GSDV-2435391920 

03:32:43.53 

-27:49:19.20 

4.763 

-0.6 

± 

0.1 

GSDV-2401153550 

03:32:40.11 

-27:45:35.50 

4.773 

1 

0 

CO 

± 

0.1 

GSDV-2219353310 

03:32:21.93 

-27:45:33.10 

4.788 

-0.2 

± 

0.1 

ERSV-2052630041 

03:32:05.26 

-27:43:00.41 

4.804 

-0.2 

± 

0.1 

GSDV-2426693897 

03:32:42.66 

-27:49:38.97 

4.831 

-0.2 

± 

0.2 

S33166 

03:32:58.38 

-27:53:39.58 

4.40 

-0.3 

± 

0.2 

N12138 

12:36:42.25 

62:15:23.25 

4.414 

-0.2 

± 

0.1 

N23791 

12:37:20.58 

62:11:06.11 

4.421 

1 

0 

± 

0.1 

S31908 

03:32:54.04 

-27:50:00.81 

4.43 

-0.4 

± 

0.1 

N13279 

12:36:46.16 

62:07:01.83 

4.444 

-0.4 

± 

0.1 

N24628 

12:37:23.57 

62:20:38.72 

4.502 

-0.6 

± 

0.2 

N28987 

12:37:19.69 

62:15:42.46 

4.53 

-0.1 

± 

0.1 

N12849 

12:36:44.68 

62:11:50.62 

4.580 

-0.2 

± 

0.1 

N31130 

12:37:57.51 

62:17:19.10 

4.680 

-0.4 

± 

0.1 

S6294 

03:32:14.50 

-27:49:32.69 

4.74 

-0.2 

± 

0.1 

S32900 

03:32:57.17 

-27:51:45.01 

4.76 

-0.5 

± 

0.1 

S1745 

03:32:05.26 

-27:43:00.42 

4.80 

-0.1 

± 

0.1 

S3792 

03:32:10.03 

-27:41:32.65 

4.81 

1 

p 

CO 

± 

0.1 

S1669 

03:32:05.08 

-27:46:56.52 

4.82 

-0.2 

± 

0.1 

N23039 

12:37:18.07 

62:16:41.72 

4.822 

1 

0 

CO 

± 

0.3 

N6738 

12:36:23.56 

62:15:20.30 

4.889 

-0.2 

± 

0.2 

N6333 

12:36:21.94 

62:15:17.12 

4.890 

-0.2 

± 

0.1 

S20041 

03:32:33.48 

-27:50:30.00 

4.90 

-0.1 

± 

0.1 

S23745 

03:32:44.07 

-27:42:27.43 

4.923 

-0.5 

± 

0.1 
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